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. The production of the initial state of the QGP in very high-energy AA colhsions is discussed 

CN , within the framework of perturbative QCD and saturation. The next-to-leading order compu- 
J> I tation of the transverse energy of minijets is reviewed. Saturation of parton production, conjec- 
tured to occur at a dynamically determinable perturbative scale, leads to estimates of the initial 
densities. The final state multiplicities are predicted by assuming an isentropic hydrodynamical 
further evolution. Comparison with RHIC data is shown. 
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^ , 1. INTRODUCTION 
> . 

^ i Particle production at the early stages, r ~ l/pr, in ultrarelativistic heavy ion colli- 
CN I sions is dominated by processes of large transverse momentum scales pt- The partonic 
multiplicities Naa{pt > Po, V^, ^v) and transverse energies E^^{pq, y/s, Ay) produced 
in processes above some scale po 3> Aqcd into a rapidity window Ay are calculable 
in perturbative QCD (pQCD) through the formalism of collinear factorization. Subse- 
^ ', quently, the number and energy densities due to the pQCD quanta can be estimated. 

Towards smaller transverse momenta, i.e. later in formation time, the number of pro- 
duced gluons is expected to grow to the extent that non-linearities start to dominate: a 
• new dynamically generated scale appears, further growth in the number of gluons becomes 
inhibited, and gluon production "saturates" P|J^-P^. For sufficiently high cms-energies 
k> ! and large nuclei, the saturation scale becomes perturbative, ~ 1...2 GeV at RHIC...LHC 
' ■ in central AA collisions of A ~ 200 [|n|,0]. 

Gluon saturation can be viewed to take place already in the initial wave functions 
of the colhding nuclei @J^, and gluon production around the dominant saturation scale 
describable in terms of classical fields |6HT0| . As non-perturbative gluon production cannot 
be computed from truly first principles, one may equally well try to estimate the bulk 
production of gluons in an AA collision by using the pQCD component alone but by 
extending the perturbative computation from large-p^ down to a dynamically determined 
saturation scale |Tl|-[T3[1. I shall focus on the latter, the pQCD+final-state saturation 
approach, in the following. The initial-state saturation and classical gluon fields are 
discussed by L. McLerran, R. Venugopalan and D. Kharzeev in these proceedings. See 



Oh! 



also ref. |14| for self-screening in parton production, and [15] for other studies of particle 



production with saturation. 



*also: Helsinki Institute of Physics, P.O.Box 64, FIN-00014 University of Helsinki, Finland 



2 



In the final-state saturation approach |]TT| the formation time of the system can be 
estimated from the saturation scale Psat as = 1/psat- This makes the computation of the 
initial number and energy densities possible: identifying the spatial and the momentum 
rapidities, and focussing on a volume element AVi ~ TcR^TiAy around rj = y = 0, we get 
the average initial densities 



„x NAAiPsi^t, Vs, Ay) E^^{p,^t,Vs,Ay) 



These densities serve as initial conditions for the further evolution of the system. I shall 
first discuss how the quantities Naa and E^^ are computed in pQCD with the conjecture 
of saturation, and then return to the multiplicity predictions for RHIC and LHC 

2. MINUETS AND THEIR TRANSVERSE ENERGY FROM PQCD 

In the leading twist approximation the semihard partonic collisions with ~ 1...2 
GeV are independent of each other. For an AA collision at an impact parameter b the 
average number of minijets produced into Ay with p^ > po and the transverse energy 
carried by them can be computed as 

Ar^A(b,po,v^,Az/) = TAAih)apQCB{N){po,y/~s,Ay,A) (2) 
E^^{h,po,V^,Ay) = TAA{h)a^QCB{ET){po,V^,Ay,A). (3) 

The nuclear overlap function TaaO^) = J c?^s T4(s)T4(b — s) takes care of the collision 
geometry. With the Woods-Saxon nuclear densities Taa{0) ~ A^/ttR^ 0. The first 
moment cr-pQCBiET) of the transverse energy distribution and o"pqcd(-^) of the number 
distribution of the minijets in Ay, are the quantities computable in pQCD. 

2.1. Leading order 

In leading order (LO), the semihard partonic collisions are 2^2 processes. Their 
inclusive cross section can be obtained through coUinear factorization, 

^^AA^kl+X da'^^''^ 

dpj^dyidy2 V ki dt 

with the momentum fractions Xi^2 = '^(^^^^ + 6^^^); the sub-cross sections 
and the factorization/renormalization scales Q ~ pt- For the 1...2 GeV minijets at mid- 
rapidity {yi or y2 ~ 0) at y/s = 200. ..5500 GeV of interest here, typically x ~ 10^^.. .10"^. 
Due to the dominance of gluons at small values of x, the gluons also clearly dominate the 
production of minijets. 

Nuclear effects in Eq. are taken into account through nuclear shadowing of parton 
distributions: fi/A{x,Q) ^ fi{x,Q). For this purpose, we apply the EKS98 nuclear mod- 



ifications |]I6| which are based on a DGLAP analysis with data from deep inelastic lA 
scattering and the Drell-Yan process and sum rules as constraints. 

The transverse energy distribution of the minijets with p^ > Po and y G Ay is [H 



da 



dEi 



dpTdyidy2S2^{pi,P2)^- — / . (5) 
PoAy J 2 dpTdyidy2 
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The key element above is the "measurement function" for in Ay, 

S2'^{Pi,P2) = 5{Et- [ei + e2]pT)0(PT > Po), (6) 

where = Q{yi G Ay). This measurement function on one hand keeps track of which 
partons fall in Ay and on the other hand ensures that the collisions are hard enough to 
allow for a perturbative treatment. The first moment of the distribution becomes then 

da 

poAy Jpo "^^^ J^y J -^y^^^ dpTdyidy 
Similarly, for the number distribution of the semihard minijets in Ay, the measurement 
function is S2 {pi,P2) = S{N — [ci + e2])0(pT > Po), and the first moment 



apQCD{ET)iPo, Vs, Ay) = [ dErEx-^ = f dpr f dyi [ dy2PT-, — , -(7) 

J dEr nn.Ai, Jm JAv J dprdyTdyo 



apQCD{N){po,y/s,Ay) = [ dNN^ = [ dpr f dyi [ dy; 

J Cll\ nn.Av Jpo J Av J 



da 



Po,Ay •'po J^y J dpTdyidy2 

2.2. Next-to-leading order 

In next-to-leading order (NLO), where a ~ a^, one must include both virtual corrections 
to 2 — >■ 2 processes and real emissions, i.e. 2^3 processes. The squared matrix elements 
for the different subprocesses were computed in 4 — 2e dimensions first by R.K. Ellis and 



Sexton The subtraction procedure of getting from the squared matrix elements in 
4 — 2e dimensions to the physical (jet) cross sections was introduced by S. Ellis, Kunszt 
and Soper [0. A detailed documentation can be found in [jl9 . 



For the 2 — >■ 3 contributions with massless partons, singularities arise when the momen- 
tum of an external leg becomes soft, or when any two external particles become collinear. 
In 4 — 2e dimensions, these singularities appear as terms ~ 1/e and Similar terms 

arise also from the 1-loop graphs (interfered with the LO 2—^2). In the sum of the 
different contributions, formed for computing a cross section of a physical observable, the 
terms singular at the limit e — >■ cancel, provided that one properly defines 

• infra-red safe generalization of the measurement function S2 to 6*3 for the 3-particle 
final state: 6*3 — > 5*2 when any final state parton becomes collinear with any other 
parton in the final or initial state, or, when any final state parton becomes soft [18|. 

• NLO PDFs to absorb one 1/e singularity in the initial state (MS scheme), 

• infra-red safe choice of the renormalization and factorization scales. 

In the computation of the minijet Et production in NLO in pp and AA collisions 



20| , |2T| , we have (see also p2[) adopted the subtraction method [T^. Although the minijet 
measurement functions differ from those for observable jets, the computational procedure 
itself is similar to that of inclusive jet production |jl9|, see [0] for details. For e.g. 



the 2 — >■ 3 terms we need to numerically evaluate 6-dimensional integrals with various 
kinematical cuts. 

The infra-red safe measurement functions 82^ for the 2^2 terms and Sf"^ for the 
2 — * 3 terms are designed to answer the following question: How much Et is carried into 
Ay by minijets produced in partonic collisions where at least an amount 2po of Et is 
produced? They are 

S2'^{Pi,P2) = KEt - [eim + e2Pr2])0(m + Pt2 > 2po) 
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5{Et - [eiPTl + ^2PT2 + e3PT3])0(PTl + PT2 + Ptz > 2po) (9) 



with t/i and pTi refering to the rapidities and transverse momenta of the final state partons. 
The semi-inclusive minijet ii^T-distributions can then be defined in NLO as 



da 



poAy 



I "^^^^^'WW^^'^^^''^'^ ^ / d[PSh^^^S^^p1,p'2,p's), (10) 



d[PS], 



where d[PS]2 {dlPS]^) is the 2(3)-particle phase space volume element. The first term 
now includes both ~ and ~ contributions. The first Et moment becomes now 



c^pQCD(^r)(Po,v^,A2/)= / dETEr—- =^{ETfryl, + cr{ETf^l^, 

Jo CLnjT poAy 



(11) 



where 



aiEr)'-' 



aiEr)'-' 



d[PS]. 
d[PSl 



da 



2-»2 



d[PS]2 
da^ 
d[PS], 



ei + £2 



PT2Q{PT2 > Po 



eiPri + e2PT2 + ^sPts 



) (12) 
6(Pti +Pt2+Pt3 > 2p()|3) 



For the renormalization and factorization scales we use Q = n{Y^iPj'i)/2 with n ~ 1. 

Figure |l| shows the first moment of the minijet i?T-distribution as a function of the 
minimum p^-scale po, computed for pp collisions at y/s = 200 GeV and y/s = 5500 GeV. 
The solid curves are the full NLO result with 2-loop a, and NLO PDFs, GRV94 (left) 



a: 



and CTEQ5 (right) |2^. The dotted curves (LO) show the lowest order results 
computed with the corresponding LO sets of PDFs and 1-loop as, and the dashed curves 
are the LO results with NLO PDFs and 2-loop as- 
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Figure 1. (Jpqcd(-E't)(P07 \/^) vs. po, from |21|. For the details, see the text. 



We see that even towards the perturbatively small scales of 1 GeV at RHIC and 2 GeV 
at the LHC, the NLO results do not grow rapidly relative to the LO (or LO') results. 
This speaks for the applicability of pQCD at the semihard scales. If a K-fa.ctoT is defined 
as NLO/LO, we notice the stability of K within po = 1...2 GeV but that quite obviously 
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the i^-factor depends on and on the PDF set used. It is also interesting to observe 
that a{Ej') at the LHC is insensitive to the choice of PDFs. 

Even though the NLO resuhs are stable relative to the ones in LO, the i^'-factors are 
rather large, and one should be concerned of the convergence of the perturbation series. 
The dependence of (t{Et) on the scale choice is very similar, quite weak, both in NLO 
and in LO The large NLO contributions mainly come from a kinematical domain 
Et < Po which is not included at all in LO. The new region is due to processes where 
two hard partons fall outside Ay and one parton with pj- < pq falls inside. To compare 
the NLO and LO results in the common kinematical region Et > po, we have excluded 
the new region Et < Po in Fig. | (left). The iT-factors now reduce dramatically: at the 
scales 1..2 GeV K ~1 relative to LO. This also speaks in favour of the applicability of 
pQCD at the semihard scales. At the same time it shows, however, that to search for a 
convergence of the perturbative series, a NNLO computation is needed. 
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Figure 2. Left: As Fig. 1 but with the region Et < Po excluded. Right: The dependence of 
a{ET) on the rapidity interval Ay. The figures are from |21]. 



The inclusive jet cross sections are known to depend on the jet cone radius R as 
dajet/dpTdy ~ CR^ + BlogR + D |T^. If such a dependence of a{ET) on Ay were 
found, the initial densities in Eqs. (P would not be well-defined but the local density 
would badly depend on Ay. Fig. || (right panel) verifies that this is not the case but that 
the NLO results for a{ET) depend in a same, linear, way on Ay as the LO results do. The 
initial energy density at r = 1/po can thus be computed as in LO before but using the 
Et obtained in NLO. Unfortunately, the NLO situation is not as straightforward for the 
number density n^. The number of partons in Ay is not an infrared-safe quantity unless 
one introduces an additional resolution scale for two nearly collinear partons. 

To conclude the pQCD part, let me note that both Naa{po) (LO) and E^'^^po) (LO 
and NLO) depend strongly on the minimum transverse momentum scale po. Therefore 
we cannot compute the initial densities and rii from pQCD alone but some additional 
(phenomenological) element of QCD must be introduced to dynamically determine the 
relevant scale Pq. In the following, saturation of produced partons [|1T| is suggested as the 
mechanism which determines the dominant scale for gluon production. 
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3. SATURATION IN PARTON PRODUCTION 

Saturation is a dynamical non-abelian feature of gluon fields. Saturation phenomena 
were first discussed by Gribov, Levin and Ryskin |^ in the context of particle production in 
pp collisions. Saturation in gluon production in AA collisions was first discussed by Blaizot 
and Mueller [|l| , then by McLerran and Venugopalan in describing the particle production 
in AA in terms of classical gluon fields @,0], also studied by Mueller |p, Kovchegov and 



Kharzeev et al |10[- These studies emphasize the saturation phenomena already in the 
initial wave functions of the colliding objects. Saturation of gluons at very small values 
of X then also saturates the actual production (liberation) of gluons. 

Saturation of gluon production can also be conjectured to arise in the following way: 
Let us consider a nuclear disc of a transverse area nR\, filled with N^{Q) gluons, all 
produced at a scale Q and correlated over some longitudinal length. The average area 
density of gluons is p{Q) ~ Nf{Q)/TiR\. The average charge within some area vri? is 

then (gs) ~ QsVN ~ Qgy/pW, and the (Coulomb-like) potential A ~ {qs)/R ~ 9s\fP- 
Saturation takes place when non-linearities in the field-strength tensor of QCD, F ~ 
dA — QsA^, become important, i.e. when dA < QsA"^. Dimensionally, dA ~ A/R ~ AQ, 
and we obtain a criterion for saturation as N^{Q) x al/Q"^ ^ 'nR\. Thus, if N^{Q) can 
be computed, the saturation scale Qsat can be obtained. 



The above is, however, at best a scaling argument. In practise, in ref. [|rT|, we have 
used a geometrical saturation criterion 

Naa{po, v^, Ay) X ^ = -rRI (14) 

where we have neglected possible effects due the running of the coupling, and have not 
included any group theoretical factors explicitly. The correlation length in rapidity has 
been taken to be Ay = 1. In other words, we extend the computation of the number 
of minijets, Naa{Po) discussed in the previous section, down to scale Pq = Psati^/s, A), 
which is determined from Eq. (|T^), and at which the produced gluons start to overlap 
transversally. The multiplicity at saturation, NAA{Psa.t), should then give a good first 
estimate of the total parton (gluon) production in central AA collisions. Notice that it is 
the integral over pt that matters and the contribution from pt < Psat to the integral is 
included through the saturation criterion. 

Figure ^ (left) shows the average number partons produced into Ay = 1 with pt > Po 
in a central AA collision oi A = 208 at y/^/A = 20 GeV (SPS), 200 GeV (RHIC) at 5500 
GeV (LHC). The curve labelled as "saturation" is plR\ from Eq. (0), and the saturation 
scale at each energy can be read off from the intersection points of the curves. 

It is quite interesting to notice that at the scaling limit (neglecting the small-x rise of the 
gluon densities, and the effects of the phase space) crpQCD(po) ~ Po^, and the saturation 
criterion (0) results in psat ~ A^^^ and in the multiplicity as A^AA(Psat) ~ A - instead of 
the A^^^ scaling typical for hard processes. 

With realistic gluon densities (GRV94LO |]23|) including the EKS98 shadowing effects 



1^, and an overall K = 2 to roughly simulate the NLO effects, we have obtained in [|TT 
the scaling laws for Psi^tiAVs), iVAyi(Psat, Ay) and E^^(psat,Ay) presented in Table 1. 



Recently, we have also shown analytically how these scaling laws arise |T3| from the initial 
gluon densities, and that A^AA(Psat) ~ ^«s(Psat)2^5'A(2psat/v^,Psat)5 where the explicit 
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Figure 3. Left: Determination of the saturation scale. Right: Charged particle multiplicity 
dNch/dy in central AA collisions for A = 208, 136, 64, 32 and 12, as a function of ^/s. See the 
text for details. The figures are from [11|. 



power of now depends on the one included explicitly in the saturation criterion (|T4t). 
At the collider energies, we see that the saturation scales are Psat ~ 1 ... 2 GeV ^ Aqcd 
for the heavy nuclei A ~ 200. 



Table 1 

The scaling laws for the saturation scale psa.t, initial multiplicity A^aa (Psat > Ay) and transverse 
energy E^^{psat, Ay), and the subsequent initial conditions for QGP at r^. 



A = 208, Ay = 1 


RHIC 


LHC 


scaling in A and ^/s 


yi/AGeV 


200 


5500 




Psat(A,v^)/GeV 


1.13 


2.13 


0.208A0-i28(^)0-i9i 


A'AA(Psat, Ay) 


1440 


5140 


1.383A°-922(yi)0.383 


£;^^(Psat,A|/)/GeV 


2360 


17000 


0.386^1043(^^^0.595 


n = 1/psat (fm) 


0.17 


0.093 




nj/fm~^ 


59.8 


401 


0.370^0.383(^)0.574 


ei/GeVfm-3 


98.2 


1330 


0.103A0-504(yi)0.786 


Ti/GeV (w. gluon d.o.f.) 


0.62 


1.19 


0.111A0-126(^)0.197 



The average initial conditions for the QGP obtained for central AA collisions with 
A = 208 for RHIC and LHC are also summarized in Table 1. The initial system is 
strongly gluon dominated, ~ 90 % of partons produced at Psat are gluons. We also note 
that from the point of view of the computed average densities ej and rii, the system looks 
thermal: the energy per gluon is as in an ideal gas, E^'^^psat) /NAA{Psat) ~ 2.7Tj, where 
Ti is the temperature of an ideal massless gluon gas at the computed energy density e^. 
Therefore, it seems that for thermalization gluon multiplication is not necessary, and early 
thermalization, perhaps already at r^, is possible (see also A. Mueller in these proceedings). 
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4. FROM INITIAL TO FINAL MULTIPLICITIES 

In between the primary production stage and decoupling, there is an expansion stage. 
A possible description of this stage is given in terms of relativistic hydrodynamics. En- 
couraged by the initial conditions for the QGP computed from the pQCD+saturation 
approach above, we assume thermalization at = 1/psat- The rapidity density of entropy 
in a longitudinally boost-invariant system is Si ~ 3.6A^^^(psat)5 counting only the gluonic 
degrees of freedom. In an isentropic expansion entropy is conserved, so Si = Sf ^ 4A^y, 
including only pions in the final state. To a first approximation, we can estimate the 
final state charged-particle multiplicity in central AA collisions directly from the initial 
conditions as 

_ 2 r |l.24A0-922(^)0.383 if g, computed from A^A^lPsat) 

iV,h = dN,i,/dy ~ 3 A^/ - j |l.l6A0-92(yi)0.40^ jf computed from E^^feat) 

This prediction is shown in Fig. ^ (right). As seen in the figure, since the system looks 
thermal from the very beginning, the scaling laws ([T5|) obtained are very close to each 
other. This illustrates the basic idea in multiplicity predictions from pQCD-fsaturation. 
For a more detailed study with y^-dependent i^'-factors, transverse profiles, transverse 
expansion effects, more detailed EoS, detailed treatment of decoupling, more complete 
list of hadrons, resonance decays, computation of particle spectra and centrality cuts, see 



refs. 25,26|, and V. Ruuskanen and U. Heinz in these proceedings. 



5. COMPARISON WITH RHIC DATA 

The above were predictions about one year before the first data from RHIC. The PHO- 
BOS experiment measured the charged-particle multiplicity in Au-Au collisions both at 
= 56 AGeV and at = 130 AGeV in Summer 2000 Fig. | from H shows the 
first PHOBOS data (6% centrality cut) by the filled circles, and the pQCD+saturation 
(EKRT) prediction [|lT| by the thick solid line. The HIJING predictions [^] with and 



without jet quenching are drawn with the solid lines. The EKRT prediction shown is 
from Eq. ([ISD with the number of participants Apart = '^A, and an approximate factor 
0.9 to account for the conversion y ^ rj at y ^ 0. 

I have also added into the figure the latest measurement by PHOBOS at y/s = 200 
AGeV (filled triangle) |3^. Considering the theoretical uncertainties, the data confirms 
the EKRT prediction amazingly well both in absolute magnitude and in the scaling with 
^/s in the RHIC energy regime. It should be emphasized that once the effective constants 
in the saturation criterion ( p!^ have been verified by comparison with the data at some 
cms-energy, the pQCD+saturation approach gives a definite prediction for the multiplic- 
ities at other energies, and also for central collisions of other nuclei. As observed in Fig. 
^, HIJING1.35 with a soft and a hard {po = 2 GeV) particle production components, 
predicts too large a multiplicity at y/s = 200 AGeV. For further discussion, see [^, and 
X.-N. Wang in these proceedings. 

Centrality dependence of the charged-particle multiplicity has been suggested as a 
further challenge for the different models ||2^. In the pQCD-|-saturation approach the 
dependence of the centrality of the AA collision can be studied by making the saturation 
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Figure 4. Left: dA'ch/'^??/(0.5A'part) as a function of y^. The prediction from pQCD+saturation 
approach |0] is shown by the thick solid line. The H1J1NG1.35 prediction is shown by the 
solid lines both for pp and for Au+Au with jet quenching (upper curve) and without (lower). 
The PHOBOS data |27| is shown by the filled circles. The figure is from ||2^, I have added the 



new PHOBOS data point at ^/s = 200 AGeV ||3^ and emphasized the saturation curve. Right: 
The same quantity as a function of the number of participants for ^/s = 130 and 5500 ylGeV as 



predicted by the pQCD+saturation approach |13,12|. The data is from PHENIX [32|. 



criterion local in the transverse plane. Noticing in Eq. (|T^) that Naa/'^R\ is the average 
transverse density of produced partons, the saturation criterion generalizes as 

cLNaa 



T^(b - s)r4(s)apQCD(iV)(Po, v^, Ay, A) = pl/n 



(16) 



data 



for an AA collision at an impact parameter b. The saturation momentum needs to be 
determined at each transverse location s, and the multiplicity of produced partons is then 
NAA{b) = J (i^spsat(b, s, ^/s, A). At sufficiently large values of s or b, psat ~ Aqcd, where 
we should not trust the pQCD calculation anymore. To avoid this, we have considered 
particle production only from the region where Psat > 0.5 GeV [|12|. Clearly, the saturation 
model [|l^] is applicable in the region where dominantly Psat ^ Aqcd- 

In Fig. ^ (right), we show the comparison of the prediction |T^ against the PHENIX 
Also a prediction for the LHC is shown. The RHIC data suggests that the 
pQCD+saturation model works for central and nearly central collisions. However, for 
non-central collisions at h>RA and A/'part^SOO, one obviously moves outside the validity 
region of the model; the dominant saturation momenta become too small, and particle 
production becomes overestimated. 

I would like to emphasize, however, that although the centrality dependence does pro- 
vide constraints for different models, the effect is still fairly modest, ~ 20 . . . 30% in the 
data at A^part ~ 100, and that in the absence of a direct measurement of Apart the systematic 
errors can be large. A more dramatic effect is predicted from the a/s- dependence of the 
charged-particle multiplicity of central collisions: as seen in Fig. | (right), A^ch/0.5Apart 
at the LHC should be almost 4 times that at RHIC [jl3l . The A-scaling of different models 



will be best tested by data from central AA collisions with different A instead of varying 
b at fixed A. 
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